Three Phases in the 3D Abelian Higgs Model with Nonlocal Gauge Interactions 
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We study the phase structure of the 3D nonlocal compact U(l) lattice gauge theory coupled with a 
Higgs field by Monte Carlo simulations. The nonlocal interactions among gauge variables are along 
the temporal direction and mimic the effect of local coupling to massless particles. In contrast to 
the 3D local abelian Higgs model having only the confinement phase, the present model exhibits the 
confinement, Higgs, and Coulomb phases separated by three second-order transition lines emanating 
from a triple point. This suggests that electron fractionalization phenomena in strongly-correlated 
electron systems may take place not only in the Coulomb phase but also in the Higgs phase. 



Some recent experiments in strongly-correlated elec- 
tron systems indicate breakdown of Landau Fermi liq- 
uid theory. Also, for certain class of quantum phase 
transitions, it is argued that Ginzburg-Landau theory 
is not applicable [1]. In these cases, the crucial point 
is the change of particle picture: quasiparticles are no 
longer electrons but certain new objects carrying frac- 
tional/exotic quantum numbers. These objects may be 
"constituents" of the original electrons. For example, 
each electron in a high-Xc cuprate may be a composite 
of a holon and a spinon. Anderson [2] argued that the 
anomalous properties of the normal state may be under- 
stood by regarding these holons and spinons as quasipar- 
ticles. This phenomenon of dissociation of electrons is 
called the charge-spin separation (CSS). Another exam- 
ple is Jain's composite-fermion picture of the fractional 
quantum Hall effect (FQHE) [3], where each electron is 
viewed as a composite of a so-called composite fermion 
and fictitous magnetic fluxes. Jain argued that when 
these fluxes cancel the external magnetic field partly, the 
FQHE is understood as the integer QHE of composite 
fcrmions in the residual field. 

To investigate these possibilities of fractionalization of 
electrons, knowledge and methods of gauge theory are 
useful. In fact, in path-integral formalism, one can intro- 
duce a U(l) gauge field as an auxiliary field that binds 
the constituents of each electron, and study the dissoci- 
ation phenomena as a confinement-deconfinement phase 
transition (CDPT) of this gauge dynamics. For the CSS, 
by starting from the t — J model and using the hopping 
expansion, we have argued that a CDPT takes place at 
certain critical temperature T = less an d the CSS is 
possible below Tqss [4]- For the FQHE, by introducing 
a fermion (we called it chargon) for a composite fermion 
and a boson (fluxon) for fluxes, we have again supported 
the CDPT [5] at T = Tpps below which chargons and 
fluxons separate. We called this phenomenon particle- 
flux separation (PFS). Furthermore, when fluxons Bose 
condense, Jain's idea of FQHE is realized. 

For the three-dimensional (3D) compact U(I) gauge 



system with local interactions without coupling to matter 
fields, Polyakov showed that the system stays always in 
the confinement phase due to the instanton condensation 
[6]. However, the issue of whether a CDPT is possible for 
a 3D (spatial 2D at T — 0)U(1) gauge system with mat- 
ter fields is still controversial [7] , mainly due to the lack 
of convincing methods for investigating a system with 
massless (gapless) matter fields whose effect is crucial for 
the phase structure. In the previous Letter [8] we intro- 
duced a 3D compact U(l) pure lattice gauge theory with 
nonlocal gauge interactions along the temporal direction 
which simulate the effect of gapless matter fields in the 
parity-symmetric phase [9]. By numerical simulations, 
we found that the CDPT occurs at a certain gauge cou- 
pling. The existence of the deconfinement phase supports 
the fractionalization phenomena like the CSS and PFS. 

To push forward this gauge theoretical approach in 
more realistic manner, one needs to include another mat- 
ter field, i.e., a Higgs field, coupled to the gauge field lo- 
cally. Actually, in the slave-boson (SB) t-J model, there 
appear fermionic spinons, bosonic holons, and a U(I) 
gauge field. In the flux state of the mean-field theory 
[10], spinons are regarded as massless Dirac fermions in 
a parity doublet. Such a spin system is studied also as an 
algebraic spin liquid [11]. To examine the effect of these 
Dirac fermions on the phase structure, etc., we perform 
path integration over the fermions formally, and approxi- 
mate the resulting fermionic determinant by the nonlocal 
gauge interactions introduced in Ref. [8] . To describe the 
bosonic holons one may introduce a U(l) Higgs field with 
the fundamental charge, its amplitude being fixed as \/5, 
where S is the holon density. In the chargon-fluxon model 
of the FQHE, integration over fermionic chargons cer- 
tainly generate nonlocal gauge interactions, and bosonic 
fluxons may be treated as a Higgs field coupled to this 
gauge field. 

In this Letter we extend the previous study [8] by cou- 
pling a U(l) Higgs field to the gauge field locally, and 
investigate the phase structure of this nonlocal Abelian 
Higgs model ( AHM) . Since the ordinary (local) 3D AHM 
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in the London limit (no radial fluctuations) has only the 
confinement phase [12,13], it is interesting to see if the 
nonlocal interactions among the gauge field are capable 
of generating a phase transition into the Higgs phase. 
In the SB t-J model and an algebraic spin liquid with 
doped holes, the possible Higgs phase implies a super- 
conducting state. In the FQHE regime, the Higgs phase 
is indispensable to achieve Jain's idea [3] since a con- 
densation of fluxons is necessary to cancel the uniform 
external magnetic field partly. 

We consider the nonlocal AHM defined on the cu- 
bic lattice of the size V = N n NiN 2 with the peri- 
odic boundary condition. The compact U(l) gauge field 
U xft = exp(i6 xft ) (— 7r < 9 xf i < 7r) is put on each link 
(x,fi) = (x,x + fi) (ji = 0,1,2), and the Higgs field 
4> x = exp(if x ) (— 7T < ip x < 7r) on each site x. The 
partition function Z is given by 
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A = A G + A H . (1) 



Aq is the gauge part of the gauge-invariant action A [14] , 
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where g is the (inverse) gauge coupling constant, V X ^ T 
is the product of U x/J , along the rectangular (x,x + i,x + 
i + t0,x + t6) of the size (1 x r) in the (i — 0) plane 
(i = 1,2), and c T is the "nonlocal coupling constant". A$ 
is the spatial single-plaquette coupling of U x ^s. Ah is 
the local Higgs action, 



Ah = {^x+^Ux^C, 



C.C. 



(3) 



where n is the Higgs coupling. In the SB t — J model, 
k cx S as explained. 

In the previous study [8] of the pure gauge system 
(k = 0), we considered (i) power-law decay c T = 1/r 
simulating the effect of massless relativistic matter fields 
and (ii) exponential decay c T = cxp(— tot) mimicking 
massive matter fields. We found a second-order CDPT 
for the case (i); the critical coupling is g c = 0.10 <~ 0.12 
for A = 1, while no CDPT for the case (ii). We also 
studied (iii) no decay case c T =const for nonrelativistic 
fcrmions in which the CDPT was observed at very small 
coupling like g c = 0.03 for V = 24 3 . Then one can nat- 
urally expect that if the Higgs transition exists in the 
power-law decay case, it also exists in the no-decay case. 
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FIG. 1. Phase diagram of the 3D nonlocal AHM. There are 
three transition lines, which are determined by the location 
of the peak of C for N = 16. 



In the following numerical studies of the gauge model 
Eq.(l), we concentrate on the power-law decay case (i) 
with A = 1 on the cubic lattice V = N 3 up to N = 
24. We first measure the "internal energy" E and the 
"specific heat" C defined as 
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^((A-{A)f). 



(4) 



In Fig.l we present the phase diagram in the g — k plane 
determined by the measurement of E and C . There are 
three transition lines separating three phases which we 
call the confinement, Higgs, and Coulomb phases. (We 
shall identify each phase later by measuring the instanton 
density and the Higgs boson mass.) 

Let us first consider the confinement-Coulomb tran- 
sition line. It includes the CDPT at n = previously 
found in Ref . [8] . Thus the deconfinement phase observed 
in Ref. [8] is identified as the Coulomb phase. We mea- 
sured E and C along a fixed k and found that E exhibits 
no hysteresis and the peak of C develops as the system 
size A^ 3 increases just as in the n — case [8]. These 
are the typical signals of a second-order phase transition. 
We note that the location of the peak of C shifts into the 
smaller g (higher-T) direction for larger TV, just the op- 
posite to the case of local interaction. It is because the 
new gauge interactions are introduced as N increases, 
and these terms favor to suppress fluctuations of U xil [8] . 

Next, let us consider the Coulomb-Higgs transition. In 
Fig. 2 we present C vs. k for g = 1.0. As N increases, 
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FIG. 2. Specific heat C vs. k for g — 1.0. The peak de- 
velpes as N increases. 
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FIG. 3. Total specific heat C (upper) and its Higgs part 
Ch (lower) vs. k for g — 0.3. 
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FIG. 5. C vs. g. For k — 1.0, C has a peak that develops 
as N increases. For k = 2.5, C has no peak. 



the peak of C also develops. Again, the location of the 
peak shifts to smaller k because of the nonlocal action 
Aq. To confirm that this transition is of second order, 
we apply the finite-size scaling hypothesis [15] to C of 
Fig. 2. We fit C as a function of N and e = (« — Koo)/^oo, 
where n x is the critical Higgs coupling at N — > oo, in 
the form of C N (e) = N^^^N^^e) with certain scaling 
function <fc(x). We get a quite good fitting for N — 8 ~ 24 
with v = 0.75, Koo = 0.23 and a = 0.12. 

In Fig. 3 we present C for a smaller value of g, g = 0.3. 
As N increases from N — 8 to 20, C gets sharper but 
the height of its peak becomes lower. To understand this 
interesting behavior, we measured the "specific heats" of 
the gauge part and Higgs part, Co — ((Ao ~ {Aq)) 2 )/V 
(O =G and H), separately [13]. Ch in Fig. 3 has a de- 
veloping peak as in the usual second-order phase transi- 
tion. On the other hand, both Cg and the cross term, 
C — Ch — Cq, are smooth and reduce as N increases. 
As explained above, the gauge specific heat Cg has the 
peak which shifts to left (smaller g) for increasing N . 
Then the tail of Cg observed at a fixed g on the right 
of the confinement- Coulomb transition line reduces as TV 
increases. This reduction gives rise to the above behavior 
of C in Fig. 3 at g — 0.3 which is close enough to the tran- 
sition line at g ~ 0.20 (See Fig.l). For larger N(> 24), 
the peak of C starts to develop as it should be since the 
Higgs part, the relevant part of the Coulomb-Higgs tran- 
sition, starts to dominate over Cq and the cross term. 
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FIG. 4. Specific heat C,C H ,C G vs. g for k = 0.25. C has 
two peaks reflecting each peak of Ch and Cg- 



Near the triple point, both Cg and Ch are relevant. 
In Fig. 4 we present C vs. g for k = 0.25. We choose 
k = 0.25 as a horizontal line because it crosses the two 
transition lines near the triple point. Each peak of Cg 
and Ch develops with N simultaneously as expected. 

Next, let us consider the confinement- Higgs transition. 
In Fig. 5 we present C vs. g for n — 1.0. At g = 0, one 
can exactly integrate over U xfl to obtain Z = [Iq(2k)] 3V , 
where Iq is the modified Bessel function. Thus Z is an 
analytic function of k, showing that there are no tran- 
sitions along the line g = 0. We measured E and C at 
several fixed k's by varying g to find that the signal of 
second-order phase transition is getting weaker as n in- 
creases. The phase transition line seems to terminate at 
the critical point, the location of which is roughly esti- 
mated as g — 0.03, k — 1.75. This critical point is viewed 
as a nonlocal version of the complementarity between the 
Higgs and confinement phases discussed in Ref. [16]. 

In order to characterize these three phases, let us cal- 
culate the instanton density p x [17]. p x is the strength 
of monopoles sitting on the site of the dual lattice x + 
(l+2+3)/2, which measures the disorderedness of gauge- 
field configurations. In Fig. 6 we show the average density 
p = (\p x \) for g — 0.1 and g — 0.5. In Fig. 6(a), p is fairly 
large at the small k region, and as crossing the transition 
point, it decreases rapidly. In Fig. 6(b), p is already very 
low even in the small k's, and it decreases further across 
the transition point. Such behavior of p and the gen- 
eral relations p (confinement) > p (Coulomb) > p (Higgs) 
lead us to identification of each phase as shown in Fig.l. 

For the AHM, one may study vortices of the Higgs field 
and their trajectories which form closed loops or termi- 
nate at (anti-)instantons [18]. The confinement phase 
is a plasma phase of instantons [6], i.e., the CDPT is 
characterized by dissociations of instanton-anti-instanton 
dipoles connected by vortex lines into a system of in- 
stanton plasma [18]. In Fig. 6 we also plot the density 
of isolated instantons Single by subtracting the nearest- 
neighbor dipole instanton pairs in the temporal direction. 
Fig. 6 clearly shows that jingle survives only in the con- 
finement phase as it should be. 
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FIG. 6. Instanton density p for N = 16 vs. k for (a) g = 0.1 
and (b) g = 0.5. The crosses are the density p s i ng ic of iso- 
lated instantons. Dashed lines indicate the transition points. 
p(g = 0.5) is smaller than p(g — 0.1) by the factor ~ 10~ 4 . 



Finally, let us consider masses of gauge-invariant oper- 
ators. In Fig. 7, we show the masses Mh (inverse correla- 
tion length) of Rc(cf> x+tl U xtl 4> x ) and Im(<^ x+jU [/ XjU </> x ) vs. 
k for g = 0.4 [19]. The two M H of Fig.7 exhibit the be- 
havior across the transition point similar to that of the 
4D local AHM [20]. We note that the 4D local AHM 
has a phase structure similar to the present model; three 
phases, a triple point and the complementarity. However 
at present, the confinement-Higgs transitions in the 4D 
AHM is believed to be of first order, which is in contrast 
to the present model. 

In conclusion, we numerically studied the 3D nonlocal 
AHM and found that the nonlocal interactions generate 
all the known phases of gauge theory. As explained, the 
existence of the Higgs phase supports both the supercon- 
ducting state of the t— J model and the FQHE as electron 
fractionalization phenomena like CSS and PFS. In fact, 
the high-T c cuprates near T = enter into the super- 
conducting state above some critical doping 5 C - This is 
interpreted in Fig.l as entering into the Higgs phase as k 
increases since n <x 5. For a parity-preserving QED3 cou- 
pled with a Higgs field, the result of the present model 
implies that a Higgs phase transition occurs for a suffi- 
ciently large number of fermion flavors since g is supposed 
to be proportional to the number of massless relativistic 
fcrmions. 
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FIG. 7. Masses Mh(N = 16) of real and imaginary parts 
of ((frx+^Uxufix) vs. k for g = 0.4. Mh of the real part has its 
minimum at the transition point as in the 4D local AHM. 
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